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ABSTRACT 


DNA barcoding of two marine tintinnids, Eutintinnus rectus and Schmidingerella arcuata, was performed using 
four samples collected from different sites in the north-eastern coast of South Korea. The loricae morphology was 
observed by light and scanning electron microscopy. Molecular data were analyzed using five nuclear ribosomal DNA 
markers (18S, ITS1,5.8S, ITS2, and 28S genes) and one mitochondrial marker (CO/ gene). The intraspecific pairwise 
differences of E. rectus and S. arcuata in the CO1 gene were 0.0-0.2% and 0.0-0.6%, respectively, while there were 


no differences in the 18S rDNA sequences. 


Keywords: CO1, DNA barcode, Eutintinnus rectus, planktonic ciliates, (DNA, Schmidingerella arcuata 


INTRODUCTION 


Tintinnids are globally distributed and play a key role in ma- 
rine ecosystems as an intermediate link in the microbial food 
webs (Stoecker and Capuzzo, 1990; Lynn, 2008). Moreover, 
they can be used as bioindicators of water quality and envi- 
ronmental pollution (Kim et al., 2012; Rakshit et al., 2017). 
One of the most commonly used markers for DNA barcoding 
of ciliates, is the 18S rDNA gene (Pawlowski et al., 2012; 
Stoeck et al., 2014). However, despite its many advantages, it 
lacks the sufficient resolution to discriminate closely related 
and sibling species (Pawlowski et al., 2012). This makes the 
use of alternative markers necessary, such as the mitochon- 
drial CO/ gene, which is capable of increasing the resolution 
at and below the species level (Jung et al., 2018; Park et al., 
2019). In this study, we investigated the morphology and phy- 
logeny of two tintinnids, Eutintinnus rectus and Schmiding- 
erella arcuata, collected from different locations of the East 
Sea/Sea of Japan. The molecular study was based on six gene 
markers, which were used to test and compare their ability 
to resolve the phylogenetic relationships within these marine 
tintinnid ciliates. The aim of this study is to increase the taxon 
sampling of the tintinnids, providing two reference sequences 


from different genera. We expect that our results will contrib- 
ute to the study and better understanding of population genet- 
ics, as well as to the differentiation of closely related species. 


RESULTS AND DISCUSSION 


Eutintinnus rectus (Fig. 1A, B) and Schmidingerella arcuata 
(Fig. 1C, D) were isolated from water samples from the East 
Sea/Sea of Japan, collected in April 2017 from four locations 
on the north-eastern coast of South Korea, i.e., Munam port 
(38°17'46.2"N, 128°33'1.1"E), Daepo port (38°10'17.2"N, 
128°36'27.8"E), Namae port (37°56'40.6"N, 128°47'13.7"E), 
and Gangneung port (37°46'18.45"N, 128°57'5.64"E). The 
water temperature and salinity ranged from 11.2-13.6°C and 
30.6-31.2 psu, respectively. The morphology of the two spe- 
cies was observed and studied under a stereomicroscope and 
a scanning electron microscope (Fig. 1); the identification of 
the species follows Yoo and Kim (1990) and Kim and Choi 
(2016). A single cell from each site was isolated for DNA ex- 
traction. PCR amplification was done using the primers New 
Euk A and LSU rev4 for rDNA sequences (Sonnenberg et al., 
2007), CiCO1 Fv2 and CiCO1 Rv2 for CO1 sequences (Park 
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DNA Barcode of Two Marine Tintinnids 


Fig. 1. Eutintinnus rectus (A, B) and Schmidingerella arcuata (C, D) from Bouin’s fixed samples (A, C) and in the scanning electron 


microscope (B, D). Scale bars: A-D=50 um. 


Table 1. Pairwise sequence distance (%) matrix among the ribosomal and mitochondrial DNA genes 


1 2 3 4 
1- Erec01? (MN185005, MN185013) = 0.2 0.2 0.2 
2- Erec02°(MN185006, MN185014) 0.0/0.0/0.0/0.0/0.2/0.0 = 0.0 0.0 
3- Erec03° (MN185007, MN185015) 0.0/0.0/0.0/0.0/0.0/0.0 0.0/0.0/0.0/0.0/0.0/0.0 = 0.0 


4- Erec04¢ (MN185008, MN185016) 0.0/0.0/0.0/0.0/0.0/0.0 


0.0/0.0/0.0/0.0/0.2/0.0 0.0/0.0/0.0/0.0/0.0/0.0 = 


5 6 7 8 
5- Sarc01? (MN185009, MN185017) = 0.2 0.2 0.4 
6- Sarc02° (MN185010, MN185018) 0.0/0.0/0.0/0.0/0.0/0.0 = 0.0 0.6 
7- Sarc03° (MN185011, MN185019) 0.0/0.0/0.0/0.0/0.0/0.0 0.0/0.0/0.0/0.0/0.0/0.0 = 0.6 


8- Sarc04¢ (MN185012, MN185020) 0.0/0.0/0.0/0.0/0.0/0.0 


0.0/0.0/0.0/0.0/0.0/0.0 0.0/0.0/0.0/0.0/0.0/0.0 = 


18S, ITS1, 5.8S, ITS2, 28S, and combination of the five genes (lower left) and mitochondrial CO1 gene (upper right) of the species Eutintinnus rectus 


(Erec01-Erec04), and Schmidingerella arcuata (Sarc01-Sarc04). 


The GenBank accession numbers for the markers are shown in parenthesis (rDNA, CO1). 


*Collected from Namae port. 
Collected from Munam port. 
“Collected from Daepo port. 
4Collected from Gangneung port. 


et al., 2019), and the internal primers 18SF790v2, 18SR300 
and 18SF1470 for 18S sequencing. The pairwise similari- 
ty matrices were constructed using 18S rDNA, ITS1, 5.8S, 
ITS2, 28S rDNA and CO/ gene sequences; they were aligned 
in Bioedit v7.2.5 (Hall, 1999) and calculated in MEGA 6.06 
(Tamura et al., 2012) with the p-distance method (Nei and 
Kumar, 2000). For phylogenetic analysis, the four COZ and 
18S rDNA gene sequences of E. rectus and S. arcuata were 
aligned with nine sequences of COJ and 51 sequences of 18S 
rDNA, each including two outgroup taxa, retrieved from the 
GenBank. Maximum likelihood and Bayesian inference anal- 
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yses were computed by IQ-tree 1.5.3 (Nguyen et al., 2015) 
and MrBayes 3.1.2 (Ronquist et al., 2012), respectively. 


Differences between ribosomal DNA sequences 

The five ribosomal DNA sequences of each species were 
completely identical, except for the 28S rDNA of E. rectus, 
which showed one nucleotide variation among the four spec- 
imens (Table 1). The four new rDNA gene sequences of S. 
arcuata isolated in the present work were identical to the 
two GenBank sequences (accession numbers JN871726 and 
JQ837815) (Agatha and Striider-Kypke, 2012; Bachvaroff et 
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Xystonella longicauda JQ408160 
Xystonella longicauda KT792933 
Xystonella acus KY290329 
Dadayiella ganymedes JX101852 
Parafavella jorgenseni MH673380 


ML/BI 


100/1.00 
100/1.00 


99/1.00] | parafavella gigantea MH673388 


Parafavella gigantea MH673409 
Parafavella gigantea MH673387 


Parafavella subrotundata MH673383 


99/1.00| Parafavella greenlandica MH673398 
Tintinnopsis fimbriata AY143560 
100/1.00 | Stenosemella ventricosa EU399538 
Stenosemella steini KT792927 
93/0.99 ı COdonellopsis americana AY143571 
-/- [L Codonellopsis morchella JQ408173 
Stenosemella nivalis FJ196074 
Codonellopsis nipponica FJ196072 
Codonella apicata EU399531 


98/0.99 
-/0.58 


100/1.00 
55/0.68 


95/1.00 
-/0.58 


MG594903 


100/0.99 


Outgroup + 
(D. appendiculata, 


95/0.98 Codonaria cistellula JQ408154 E. vannus) 
99/1.00 Codonaria cistellula JQ408167 co1 
53/0.59 Dictyocysta lepida KT792929 Be =a 


Dictyocysta elegans KY290318 
Stenosemella pacifica JN831794 
Codonellopsis gaussi JQ924053 
Laackmanniella prolongata JQ924056 
Helicostomella subulata JN831784 
Tintinnopsis parva JN831824 
Tintinnopsis tenuis JN831848 
Tintinnopsis parvula KU715771 
Tintinnopsis rapa JN831834 


100/1.00 
-/0.64 


99/0.99 


100/1.00 
99/1.00 
99/1.00 


51/0.63 


88/0.90 


100/1.00 
100/1.00 


100/1.00 
-10.57 


Schmidingerella arcuata JN871726 
Schmidingerella arcuata JQ837815 
Rhabdonella spiralis JQ408158 
Rhabdonella elegans JQ408175 
Cyttarocylis ampulla JQ408168 
Cymatocylis calyciformis JQ924046 
92/1.00 — Ptychocylis minor KY290321 
100/1.00 ; Tintinnopsis lobiancoi JN831814 
Stylicauda platensis JN831832 
Parundella aculeata KY290327 

Tintinnopsis urnula JN831851 

Tintinnopsis ventricosoides KU715776 
Tintinnopsis lacustris JQ408162 
Tintinnopsis cylindrica JN831811 
Tintinnopsis tubulosoides AF399109 
Tintinnopsis cf. cylindrica MG603624 
Eutintinnus pectinis AF399170 


83/0.82 


72/0.51 


85/0.88 


100/1.00 


Eutintinnus fraknoi JQ408157 


Amphorellopsis quinquealata JQ924058 
Strombidium cuneiforme KP260512 


Novistrombidium apsheronicum FJ876958 


Fig. 2. Maximum likelihood (ML) phylogenetic trees inferred from the 18S rDNA (A) and CO1 (B) gene sequences. The phylogenetic 
trees show the position of four East Sea/Sea of Japan isolates of Eutintinnus rectus and Schmidingerella arcuata. New sequences 
are shown in bold text. Numerical values for bootstrap support are given at the respective nodes as: posterior probability/ML and 
Bayesian inference (BI) bootstraps. Values <50% (or 0.50) are shown as ‘-’. The GenBank accession numbers of the 18S rDNA and 
CO1 gene sequences are provided after the species names. The scale bars correspond to three substitutions per 100 nucleotide po- 
sitions (18S rDNA tree) and 20 substitutions per 100 nucleotide positions (CO1 tree). 


lantic Ocean (Long Island Sound, USA) (GenBank accession 
number AF399170) (Snoeyenbos-West et al., 2002) (Fig. 2A). 
Eutintinnus rectus can be differentiated from E. pectinis in 


al., 2012) isolated from the Atlantic Ocean (coast of Mary- 
land, USA). The 18S rDNA sequences of E. rectus were 
identical to a sequence of E. pectinis isolated from the At- 
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having more number of oral teeth and stouter lorica (Yoo and 
Kim, 1990). Based on the picture of E. rectus, included within 
the supplementary material of provided by Snoeyenbos-West 
et al. (2002), the Atlantic Ocean and Korean isolates is partic- 
ularly different in the oral diameter (about 31 um vs. 42 um). 

Several studies have shown the inefficiency of the 18S 
rDNA to discriminate closely related species (Lynn and 
Striider-Kypke, 2006; Quintela-Alonso et al., 2013). A ge- 
netic variation of 0.6%, for the 28S rDNA, is considered as 
an ideal threshold to discriminate ciliate species (Stoeck et 
al., 2014). However, according to Santoferrara et al. (2015), 
neither the 28S nor 5.8S - ITS markers provide enough res- 
olution to unequivocally distinguish closely related tintinnid 
species. Another issue related to the rDNA loci was reported 
by Zhang et al. (2017), who showed that intraspecific vari- 
ations of rDNA loci (e.g., 18S, ITS1, 5.88, ITS2, and 28S) 
are non-significantly correlated when a 18S rDNA variation 
range of 0-1% was compared to the other loci in tintinnids. 
The use of multi-gene barcodes stands as a valuable tool to 
differentiate species, which has been successfully tested in 
tintinnids of the common genus Helicostomella (Santoferrara 
et al., 2015). 


Differences between CO1 sequences 

The intraspecific genetic differences in the CO1 gene of E. 
rectus and S. arcuata were 0-0.2% (0 or 1 nucleotide differ- 
ence) and 0-0.6% (0 to 3 nt), respectively (Table 1). Unfor- 
tunately, the dataset of CO7 gene sequences of marine tintin- 
nids is still limited to a few species of the genus Parafavella 
and to an isolate of Tintinnopsis cf. cylindrica (see Jung et al., 
2018; Park et al., 2019). Thus, further studies are necessary to 
confirm the utility of the CO7 barcoding for tintinnids. 

The existence of cryptic species within the tintinnids and 
their highly variable lorica morphology, are well-known is- 
sues for the identification and differentiation of species (Kim 
et al., 2010; Jung et al., 2018). Therefore, identifying tintinnid 
isolates based only on morphological criteria, becomes a dif- 
ficult task. Further studies are necessary to increase the taxon 
sampling, providing new reference sequences of the CO] 
gene, as well as rDNA loci, as suggested by Santoferrara et 
al. (2015). Only then, it will be possible to improve the reso- 
lution of morphological and molecular analyses, setting the 
boundaries of the interspecific variability, and enhancing the 
unambiguous and accurate differentiation of tintinnid mor- 
photypes. 
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